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На основе разработанной ранее модели выпадения примеси из облака в устойчиво стратифицированной атмосфере проведен 
расчет уровня загрязнения поверхности земли химическими веществами при испытании ракетных двигателей на твердом топли- 
ве и радионуклидами при взрыве ядерных зарядов. Проведено качественное сравнение теоретических прогнозов сданными на- 
блюдения. Предложены рекомендации по дальнейшему совершенствованию методики расчета переноса примесей в устойчиво 
стратифицированной атмосфере. 


Атмосфера Земли представляет собой с термо- 
динамической точки зрения чрезвычайно сложную 
систему. В масштабе планеты она сосредоточена в 
относительно малой области. Еще значительно ме- 
ньшие пространства занимают в ней живые орга- 
низмы. Но процессы, протекающие в атмосфере 
настолько сложны, что только лишь математиче- 
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ское описание тепловых и гидродинамических яв- 
лений в сильно упрощенном представлении ее сло- 
ем жидкости над поверхностью земли встречается с 
огромными трудностями. Их преодоление при ны- 
нешнем развитии вычислительной техники прак- 
тически невозможно. Поэтому приходится прибе- 
гать к упрощенным модельным представлениям. 



Технические науки 


Наиболее распространенные здесь подходы опира- 
ются на диффузионные и гравитационные меха- 
низмы осаждения [1]. Несколько в стороне, ввиду 
относительно редкой возможности проявления, 
находится механизм вымывания, согласно которой 
частицы примеси могут захватываться естествен- 
ными атмосферными осадками [2]. 

Как правило, в выброшенном в атмосферу 
облаке содержатся твердые или жидкие частицы 
вредных веществ, размеры которых меняются от 
нескольких единиц до нескольких десятков мкм. 
Важной величиной, по которой можно судить об 
уровне загрязнения поверхности земли, является 
количество о(х,у) выпавшего осадка на единицу 
площади поверхности земли. 

В работе [3] предложена модель осаждения мел- 
кодисперсной примеси на поверхность земли из 
облака, находящегося в устойчиво стратифициро- 
ванной атмосфере. Такое состояние атмосферы в 
приповерхностном слое воздушной массы может 
реализоваться при отсутствии прямых солнечных 
лучей на поверхность земли, т.е. в относительно 
тихую пасмурную погоду, вечернее или темное вре- 
мя суток. При наличии ветра необходимо еще, что- 
бы поверхность земли была гладкой. Точнее разме- 
ры неровностей на ней должны быть много меньше 
толщины слоя стратификации. 


Эволюция выброшенного облака загрязняющих 
веществ в устойчиво стратифицированной атмосфе- 
ре подчиняется законам, отличным от тех, что пред- 
сказываются теорией турбулентного переноса. В 
частности облако, сплющиваясь под действием сил 
интрузии в подобную диску фигуру, практически 
сохраняет свой первоначальный объем. При этом 
течение газа в нем турбулентное, а за пределами 
облака турбулентность весьма низка [4, 5]. Упро- 
щенные гидродинамические уравнения, описываю- 
щие течение жидкости внутри облака, имеют вид [3] 
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где к - полутолщина облака; г - радиальная коорди- 
ната; 1 - время; ѵ - единственная ненулевая компо- 
нента вектора ѵ скорости газа в облаке; г, - турбулент- 
ное напряжение; со ь - частота Брента-Вяйсаля, для ат- 
мосферы Земли й) 4 »0,01 с -1 . Турбулентное напряже- 
ние дается выражением т,=рк 2 ѵ 2 . Эмпирическая кон- 
станта к^0,4. Последнее соотношение в системе (1) 
является условием сохранения объема Коблака. 

Как показано экспериментальными наблюде- 
ниями [4] и анализом уравнений (1), облако в 
устойчиво стратифицированной атмосфере быстро 
за время порядка 1 /со ь выходит на продолжитель- 
ный (до нескольких сут.) автомодельный режим 
эволюции, при котором [3] 
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где і х - условное время начала автомодельного ре- 
жима. Числовой параметр <*., определяющий раз- 
мер (радиус) облака в плане, равен 2,63. 

Идея, принятая за основу модели осаждения 
примеси из облака на поверхность земли проста: по 
мере движения облака из него по направлению к зе- 
мле вылетают со скоростью ѵ,е частицы, падающие 
далее свободно в поле силы тяжести. Первый сом- 
ножитель ѵ. - это скорость свободного падения ча- 
стиц в ламинарной среде, ее можно рассчитать по 
формулам Стокса [6], Рыбчинского-Адамара [7] и 
др. известным выражениям. Второй сомножитель г 
имеет смысл интенсивности турбулентности и учи- 
тывает в простейшей форме хаотическое турбулент- 
ное движение частиц в облаке [3]. Значение е для 
чисел Рейнольдса Ке=10 6 ...Ю 7 составляет 0, 03. ..0, 05. 

Изменение а в точке (х,у) в момент времени ( 
можно задать уравнением [3] 

= Ѳ №М 1 )- ^( Х ~ ѵ х 1 ) 2 + (У- V ) 2 ]’ ( 2 ) 

а I 2Ѵ ѵ 


где ѵ„ ѵ у - компоненты вектора скорости ветра на 
высоте зависания облака в декартовой системе ко- 
ординат, связанной с поверхностью земли; 0[Д - 
единичная функция Хевисайда, если ^>0, 

0[Д=О, если ?<0. 

Функция Хевисайда определяет границу облака 
в проекции на поверхность земли. 

Ур. (2) записано при допущении равномерного 
распределения частиц в облаке. Это положение гру- 
бое, но оно сильно упрощает модель, так как отпа- 
дает необходимость в рассмотрении динамики ча- 
стиц в облаке. Для массы М справедливо уравнение 
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из него интегрированием при начальном условии 
М({=О)=М 0 получим 
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Таким образом, для нахождения поверхностной 
плотности а (1,х,у) имеем одно обыкновенное диф- 
ференциальное уравнение первого порядка 


<7 а 
сі ( 


ѵ,еМ 
2Ѵ 1 


%.Г 0 (0- у](х-ѵ х і) 2 + ( у-ѵ у і ) 2 


с начальным условием сх(/=0)=0. Здесь М опреде- 
ляется согласно формуле (3). 
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Поставленная задача численно решалась мето- 
дом Рунге-Кутта. Функция Хевисайда аппрокси- 
мировалась выражением 0[г] = [1+ехр(— ;де)]“‘, где 
\/Х - характерная длина, которую необходимо 
брать достаточно малой для четкого выделения 
границы облака. Согласно известной формуле 
Стокса [6], характерная скорость выпадения ча- 
стиц в ламинарной среде ѵ,=(2/9) р р <ір§/р {/и - дина- 
мическая вязкость газа облака). 

В качестве примера на рис. 1 приведен результат 
расчета поверхностной концентрации а, кг/м 2 , хи- 
мических загрязняющих веществ в виде А1 2 0,, НС1, 
Н,$0 4 и др. Численные значения физических пара- 
метров, отвечающих наземному испытанию пер- 
вой ступени межконтинентальной баллистической 
ракеты [3], взяты в следующем виде: р = 3000 кг/м 3 , 
^,=3- 10 6 м, 5=9,81 м/с 2 , ѵ=5 м/с, ѵ,=0 м/с, 
ц=1,8-10- 5 кг/ (м*с), К=16,010 6 м 3 , к=0,3, е=0,03, 
Л/ 0 = 5 • 1 0 4 кг, 1/*=10м. 



X , км 

Рис. 1 . Изолинии уровня поверхностной плотности частиц 
ег 1СГ 4 (х,у) кг/м 2 . Номера соответствия областей, огра- 
ниченных изолиниями со значениями ст, следующее: 
1) 0,87..1,0 ; 2) 0, 74... 0, 87; 3) 0,61.0,74; 4) 0,48...0,61; 
5) 0,35-0,48; 6) 0,22-0,35; 7) 0,09-0,22; 8) 0-0,09 

Расчеты также показывают, что наиболее силь- 
но влияющими на Ых,у) размерными параметрами 
являются ѵ. и ѵ„ ѵ у (точнее, безразмерные комплек- 
сы ѵ./ѵ„ ѵ./ѵ,). 

Сделаем несколько замечаний. В практическом 
приложении изложенной модели необходимо 
иметь в виду, что частицы, покинув облако, прео- 
долевают еще расстояние от облака до поверхности 
земли. Поэтому, если за точку выброса считать на- 
чало координат, результаты расчетов а по выше- 
приведенным уравнениям должны быть смещены 
на некоторое расстояние (в самом уравнении для 
а этого можно достичь, заменив х-^х+хц). Расстоя- 
ние х 0 должно определяться с учетом неравномер- 
ности скорости ветра с высотой. 

Апробация любой модели на эксперименталь- 
ных данных по радиоактивному загрязнению по- 
верхности земли в результате испытаний ядерного 
оружия или техногенных катастроф представляет 
собой достаточно сложную задачу. Сложность 
сравнения теоретических и экспериментальных ре- 
зультатов объясняется следующими существенны- 
ми факторами: 


- обилием в продуктах взрыва ядерного заряда 
радионуклидов, постоянные распада которых 
меняются в широком диапазоне; 

- большим различием в количествах образую- 
щихся при взрыве радиоактивных элементов; 

- расположением эпицентра ядерного взрыва от- 
носительно поверхности земли и вида этой по- 
верхности (тип водного пространства, характер 
грунта и пр.); 

- погодными условиями; 

- разновидностью и мощностью ядерного заряда; 

- временем проведения измерений уровня радио- 
активного заражения и др. 

По этим причинам при разработке моделей для 
расчета загрязнения окружающей среды радиону- 
клидами на первом этапе целесообразно пользо- 
ваться параметрами, имеющими смысл усреднен- 
ных величин по количеству возможных радиоак- 
тивных изотопов и интенсивности излучения по 
основным его видам: а-частицам, /3- и у-лучам. За- 
частую такой подход оправдан тем, что результаты 
измерения отражают интегральные и усредненные 
характеристики уровня радиации. 

Для ясного представления масштабности и слож- 
ности расчета заражения от радиоактивного облака, 
возникающего при взрыве ядерных зарядов, рассмо- 
трим в качестве примера испытание «Сайіе Вгаѵо», 
произведенного 1 марта 1954 г. на атолле Бикини. Кра- 
тер от взрыва получился 2 км в диаметре и глубиной 
75 м [8]. Объем выброшенного грунта составил поряд- 
ка 10* м 3 с общей массой около 150 млн т. За 3 с огнен- 
ный шар вырос до 5,5 км в диаметре. В течение мину- 
ты грибовидное облако достигло высоты 15 км, под- 
нявшись до 30 км двумя минутами позднее. За 6 мин 
облако взлетело до 40 км, через 8 мин после взрыва 
гриб достиг максимального размера в 20 км в диаме- 
тре, со стволом толщиной в 7 км, нижний его край на- 
чинался на высоте 16,5 км. «Сайіе Вгаѵо» вызвал кру- 
пнейшее в США радиоактивное заражение прилегаю- 
щих территорий и облучение местных жителей. 

В представленном ниже рис. 2 (внизу) даны из- 
меренные суммарные дозы облучения [8]. Здесь же 
приводятся результаты расчета (вверху) поверх- 
ностной концентрации загрязняющих веществ, 
выпавших на землю. Изолинии прочерчены по 
1,33-кратному падению поверхностной концентра- 
ции а от предыдущего значения. Области с кон- 
центрацией частиц ниже десятикратного падения 
от максимального значения в центральной части 
области загрязнения не приводятся. Здесь возмож- 
ность сравнения расчетных данных с результатами 
измерения сильно ограничивается отсутствием ме- 
теорологических данных на момент испытания. 
Также заведомо неизвестно количество выброшен- 
ного в атмосферу радиоактивного вещества и ха- 
рактеристики его радиационного излучения, то 
сравнение проводится только по геометрическим 
размерам области выпадения радиоактивных ве- 
ществ и темпу падения уровня радиации как функ- 
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ция от времени и расстояния от эпицентра взрыва. 
Этот шаг корректный, т.к. суммарная доза облуче- 
ния напрямую зависит от количества имеющегося 
источника радиации. 

Действительно, если А г (с -1 ) - интенсивность 
испускания у-квантов со средней (по всем возмож- 
ным источникам данного вида излучения) энерги- 
ей Е у , то количество энергии д,, испускаемое в еди- 
ницу времени числом ядер Ы у загрязняющих ве- 
ществ, равно [9] 


В формуле 


«г=\ Е Л- 


N = —М А , 

ѵ А ? 


выражающей число ядер у-излучателей А г через их 
полную М у и молярную ц у массы, число Авогадро 
И А , заменим М г на поверхностную плотность а у 
примеси. Тогда д 7 будет иметь смысл поверхност- 
ной плотности энергии у-квантов, излучаемой в 
единицу времени от выпавших на поверхность зе- 
мли загрязняющих радиоактивных веществ: 

? г =даЛ[ Дж/мЧ]. (4) 

При данном коэффициенте поглощения телом 
у-квантов к, можно найти суммарную дозу облучения 
(Дж/кг = Гр) путем интегрирования по времени и его 
поверхности выражения к,д ѵ и поделив полученный 
результат на массу рассматриваемого тела [2, 10, 11]. 
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Рис. 2. Изолинии ( внизу) суммарной дозы облучения ( рад) при 
непрерывной экспозиции продолжительностью до 96 ч 
после наземного испытания «СазІІе Вгаѵо». Вверху ~ ре- 
зультаты расчета качественного поведения поверхност- 
ной концентрации загрязняющих частиц по модели на- 
стоящей работы. Номера соответствия областей, ограни- 
ченных изолиниями со значениями о1(Т 5 кг/м 1 , следую- 
щее: 1) 1,75..2,00; 2) 1,50...1,75; 3) 1,25...1,50; 4) 1,00.2, 25; 
5) 0,75...1,00; 6) 0,50...0,75; 7) 0,25...0,50; 8) 0...0.25 

Аналогичным образом можно получить выра- 
жения для поверхностных плотностей д а , д р энер- 
гии а- и Д-излучения 


Ча = К Е а**А — > Чр =\Е р Ы А ^ (5) 

Ча Чр 

с константами и энергиями распада А я , Хр, Е а , Ер, 
поверхностными плотностями осадков сг„, а р , и мо- 
лекулярными массами <т,„ и соответствующие 
дозы облучения. 

При этом общее количество радиоактивных ос- 
адков а связано с осадками по отдельным видам 
излучения простым соотношением 

17 _ Т* , , °у 

(ч) Ча Ч Р Чу ’ 

где (ц) - средняя молярная масса радиоактивных 
осадков. 

Приведем грубую теоретическую оценку погло- 
щенной человеком дозы радиации в энергетическом 
эквиваленте в районе радиоактивного загрязнения 
атолла Бикини. Опуская индексы, указывающие на 
тип излучения, примем по порядку величин эффек- 
тивный период полураспада радионуклидов при- 
мерно равным одному году, т.е. А=0,ЗЗГ(Г 7 с -1 ; сред- 
нюю энергию излучения 2?=0,01 МэВ=1,6Т(Т 15 Дж; 
поверхностную плотность загрязняющих веществ 
о=10 мг/м 2 =10~ 5 кг/м 2 ; среднюю молярную массу 
уП=0,032 кг/моль. Тогда имеем <7=10~ 2 Дж/(м 2 -с). За 
100 ч человек массой 60 кг с площадью поверхности 
тела ~1 м 2 и коэффициенте поглощения 0,5 получит 
дозу ~30 Дж/кг=3000 рад, что находится в хорошем 
согласии с данными наблюдения (рис. 2). Взятая 
здесь величина о=10~ 5 кг/м 2 соответствует начально- 
му 4/ 0 =ЮО т радиоактивного вещества в облаке. 

В формулы (4, 5) не входят энергетические пото- 
ки от газообразных радиоактивных веществ. Для 
них расчет должен вестись по их массовой концен- 
трации У [кг/м 3 ] содержания в воздухе. Искомые 
энергетические потоки различных видов излучения 
получаются из вышеприведенных формул путем за- 
мен сг »Г„, Ор -» Ур , <т, — э К, При этом необходимо 
учитывать еще поглощающую способность самой 
атмосферы, которая оказывает довольно существен- 
ное ослабляющее действие на а- и Д-излучения. 

Все приведенные рассуждения касаются случа- 
ев, реализующихся при асимптотически больших 
временах, когда есть смысл разграничивать радио- 
активные вещества на твердые (или жидкие) вы- 
павшие осадки и газовые компоненты. На стадии 
выпадения для осадков, когда часть их остается 
еще во взвешенном состоянии в воздухе, расчет до- 
зы облучения ведется одновременно как через мас- 
совую концентрацию У мелкодисперсных частиц, 
так и их поверхностную концентрацию. 

Скорость ветра, судя по эффективному распро- 
странению облака, менялась относительно незна- 
чительно во времени и более существенно по про- 
странству. Грубый расчет дает оценку для модуля 
вектора скорости около 10 м/с. Поэтому принято 
ѵ = 1 0 м/с. Для поперечной компоненты ѵ у вектора 
скорости ветра использовались различные формы. 
Результаты на рис. 2 проведены с формой ѵ у = 1 0 _4 -х, 
где расстояние х измеряется в км. 
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Расчет с квадратичным представлением 
ѵ у =1,0— 0,01х+2,510' ~ 5 -х* несколько лучше воспро- 
изводит общую картину загрязнения радиоактив- 
ными веществами (рис. 3), особенно в районе атол- 
лов Ронгелап, Алингинае и Ронгерик, но хуже на 
больших расстояниях от эпицентра. 



X , км 

Рис. 3 . Испытание «СазІІе Вгаѵо». Картина уровня загрязне- 
ния поверхности земли при ветре с квадратичным 
профилем скорости ветра от координаты х. Соответ- 
ствие нумерации областей со значениями поверх- 
ностной плотности <у такое же, как и на рис. 2 

В проведенных расчетах моделирования рас- 
пространения радиоактивных веществ при испыта- 
нии «СазДе Вгаѵо» эффективная скорость осажде- 
ния частиц ѵ,е=1(Р 4 м/с. 
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том явления турбулентной диффузии можно зада- 
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ной ширины следа к его характерной длине соста- 
вляет около ОД. 

Выводы 

Подводя итог к проделанной работе, можно за- 
ключить, что предложенная модель может приме- 
няться к оценке уровня загрязнения атмосферы хи- 
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устойчиво стратифицированном состоянии атмо- 
сферы. 
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